AD-A151  572 


ABSTRACT 


■>  Precipitation  data  are  obtained  from  the  GWE  (Global  Weather 
Experiment)  FGGE  (First  GARP  (Global  Atmospheric  Research  Project) 
Global  Experiment)  data  sets  observed  from  December  1978  through 
November  1979.  The  FGGE  precipitation  data  are  compared  to  clima¬ 
tological  data  and  other  fields  compiled  or  calculated  from  the  FGGE 
data.  These  other  fields  are  the  U  and  V  wind  components,  heating 
fields,  velocity  potentials,  and  stream  functions.  The  precipitation 
data  are  edited  and  put  into  a  useful  global  format.  The  1979  data 
are  compared  to  climatology.  The  precipitation  fields  generally 
agree  with  more  extensive  climatologies,  though  some  important 
regional  differences  are  noted.  Large  scale  overturning  circulations 
are  colocated  with  extensive  areas  of  precipitation  in  the  tropics. 

The  analysis  of  rainfall  is  adversely  affected  by  the  sparse 
data  density  over  the  oceans  and  in  remote  areas.  In  the  future, 
the  precipitation  data  should  be  combined  with  and  scrutinized 


against  satellite  data,  especially  over  the  oceans. 


Me cession  For  mmm 
HTtis  gram 

I  DT IC  TAB  J5 

Unannounced  p. 
Justification. _ _ 

By — . 

_Distributi^7~~  ' 

— a  bl  1  i  ty~Coc3e  s 

and/or  " 
Dlst,  Special 


AFIT  RESEARCH  ASSESSMENT 


The  purpose  of  this  questionnaire  is  to  ascertain  the  value  and/or  contribution  of  research 
accomplished  by  students  or  faculty  of  the  Air  Force  Institute  of  Technology  (AU).  It  would  be 
greatly  appreciated  if  you  would  complete  the  following  questionnaire  and  return  it  to: 

AF IT/NR 

Wright-Patterson  AFB  OH  45433 

RESEARCH  TITLE-  Characteristics  Of  Precipitation  Regimes  During  First  Global  Atmospheric 

_ Research  Project  Global  Experiments _ 

AUTHOR:  _ Glenn  David  Ahrens _ 

RESEARCH  ASSESSMENT  QUESTIONS: 

1.  Oid  this  research  contribute  to  a  current  Air  Force  project? 

(  )  a.  YES  (  )  b.  NO 

2.  Do  you  believe  this  research  topic  is  significant  enough  that  it  would  have  been  researched 
(or  contracted)  by  your  organization  or  another  agency  if  AFIT  had  not? 

(  )  a.  YES  (  )  b.  NO 

3.  The  benefits  of  AFIT  research  can  often  be  expressed  by  the  equivalent  value  that  your 
agency  achieved/received  by  virtue  of  AFIT  performing  the  research.  Can  you  estimate  what  this 
research  would  have  cost  if  it  had  been  accomplished  under  contract  or  if  it  had  been  done  in-house 
in  terms  of  manpower  and/or  dollars? 

(  )  a.  MAN-YEARS  _  (  )  b.  $ _ 

4.  Often  it  is  not  possible  to  attach  equivalent  dollar  values  to  research,  although  the 
results  of  the  research  may,  in  fact,  be  important.  Whether  or  not  you  were  able  to  establish  an 
equivalent  value  for  this  research  (3.  above),  what  is  your  estimate  of  its  significance? 

(  )  a.  HIGHLY  (  )  b.  SIGNIFICANT  (  )  c.  SLIGHTLY  (  )  d.  OF  NO 

SIGNIFICANT  SIGNIFICANT  SIGNIFICANCE 

5.  AFIT  welcomes  any  further  comments  you  may  have  on  the  above  questions,  or  any  additional 
details  concerning  the  current  application,  future  potential,  or  other  value  of  this  research. 

Please  use  the  bottom  part  of  this  questionnaire  for  your  statement(s). 


CHARACTERISTICS  OF  PRECIPITATION  REGIMES 


DURING  FIRST  GLOBAL  ATMOSPHERIC  RESEARCH 
PROJECT  GLOBAL  EXPERIMENTS 

by 

Glenn  David  Ahrens 
Captain,  Air  Force,  1984,  73  pages 

Master  of  Science 
The  University  of  Utah 


CHARACTERISTICS  OF  PRECIPITATION  REGIMES 


DURING  FIRST  GLOBAL  ATMOSPHERIC  RESEARCH 
PROJECT  GLOBAL  EXPERIMENTS 

by 

Glenn  David  Ahrens 
Captain,  Air  Force,  1984,  73  pages 

Master  of  Science 
The  University  of  Utah 


ABSTRACT 


Precipitation  data  are  obtained  from  the  GWE  (Global  Weather 
Experiment)  FGGE  (First  GARP  (Global  Atmospheric  Research  Project) 
Global  Experiment)  data  sets  observed  from  December  1978  through 
November  1979.  The  FGGE  precipitation  data  are  compared  to  clima¬ 
tological  data  and  other  fields  compiled  or  calculated  from  the  FGGE 
data.  These  other  fields  are  the  U  and  V  wind  components,  heating 
fields,  velocity  potentials,  and  stream  functions.  The  precipitation 
data  are  edited  and  put  into  a  useful  global  format.  The  1979  data 
are  compared  to  climatology.  The  precipitation  fields  generally 
agree  with  more  extensive  climatologies,  though  some  important 
regional  differences  are  noted.  Large  scale  overturning  circulations 
are  colocated  with  extensive  areas  of  precipitation  in  the  tropics. 

The  analysis  of  rainfall  is  adversely  affected  by  the  sparse 
data  density  over  the  oceans  and  in  remote  areas.  In  the  future, 
the  precipitation  data  should  be  combined  with  and  scrutinized 
against  satellite  data,  especially  over  the  oceans. 


CHARACTERISTICS  OF  PRECIPITATION  REGIMES 


DURING  FIRST  GLOBAL  ATMOSPHERIC  RESEARCH 
PROJECT  GLOBAL  EXPERIMENTS 


by 

Glenn  David  Ahrens 


A  thesis  submitted  to  the  faculty  of 
The  University  of  Utah 

in  partial  fulfillment  of  the  requirements  for  the  degree  of 


Master  of  Science 


Department  of  Meteorology 
The  University  of  Utah 
December  1984 


I  H I  IMMUNITY  Ol  1  I  AH  GU  M  >1  A  1  !  NCHOOl 


FINAL  R  LADING  APPROVAL 


To  the  Graduate  Council  of  The  University  of  Utah 


Glenn  D .  Ahrens 


1  ha\e  read  the  dissertation  of  oiemi  u‘  "  "  r  K  "  * _ in  its 

final  form  and  have  found  that  (1)  its  format,  citaiions.  and  bibliographic  style  are 
consistent  and  acceptable.  (2)  its  illustrative  materials  including  figures,  tables,  and 
charts  arc  in  place:  and  (?)  the  final  manuscript  is  satisfactory  to  the  Supervisory 
Committee  and  is  reads  for  submission  to  the  Graduate  School 


Member.  Supervisory  C  ommitlet 


Approved  fot  the  Major  Department 


J  .  E  .  Ge  i  s  1  e  r 


C  hairman  t  Dear 


Approved  for  the  Graduate  Council 


James  L.  Clayton 

Dean  of  1  he  (iraduatr  Si  hot*! 


THE  UNIVERSITY  OF  UTAH  GRADUATE  SCHOOL 


SUPERVISORY  COMMITTEE  APPROVAL 


of  a  thesis  submitted  bv 

Glenn  D.  Ahrens 


I  his  t  nests  has  been  read  b\  each  member  of  the  following  supers  isorv  committee  and  bs  majont 
sote  has  been  tound  to  be  sattsfactors 


ABSTRACT 


Precipitation  data  are  obtained  from  the  GWE  (Global  Weather 
Experiment)  FGGE  (First  GARP  (Global  Atmospheric  Research  Project) 
Global  Experiment)  data  sets  observed  from  December  1978  through 
November  1979.  The  FGGE  precipitation  data  are  compared  to  clima¬ 
tological  data  and  other  fields  compiled  or  calculated  from  the  FGGE 
data.  These  other  fields  are  the  U  and  V  wind  components,  heating 
fields,  velocity  potentials,  and  stream  functions.  The  precipitation 
data  are  edited  and  put  into  a  useful  global  format.  The  1979  data 
are  compared  to  climatology.  The  precipitation  fields  generally 
agree  with  more  extensive  climatologies,  though  some  important 
regional  differences  are  noted.  Large  scale  overturning  circulations 
are  colocated  with  extensive  areas  of  precipitation  in  the  tropics. 

The  analysis  of  rainfall  is  adversely  affected  by  the  sparse 
data  density  over  the  oceans  and  in  remote  areas.  In  the  future, 
the  precipitation  data  should  be  combined  with  and  scrutinized 
against  satellite  data,  especially  over  the  oceans. 


My  thesis  is  dedicated  to  my  wife  and  family  for  the  long 
hours,  support,  and  assistance  during  the  writing  of  this  thesis. 


TABLE  OF  CONTENTS 


ABSTRACT  . 

ACKNOWLEDGEMENTS  . 

Chapter 

1.  INTRODUCTION . 

1.1  Motivation  . 

1.2  Background  . 

2.  PROCEDURES . 

2.1  Data  Acquisition  . 

2.1.1  Precipitation  Analyses. 

2.1.2  Heating  Rate  Fields  .  . 

2.1.3  Circulation  Diagnostics 

2.2  Analysis  Methods  . 

2.2.1  Precipitation  Analysis. 

2.2.2  Heating  Rate  Analysis  . 

2.2.3  Circulation  Analysis.  . 

3.  RESULTS  . 

3.1  Introduction  . 

3.2  Seasonal  Averages . 

3.3  Weekly  Averages . 

3. 4  Dai ly  Variations  . 

3.5  Hourly  Averages . 

4.  CONCLUSIONS  . 

4.1  Data  . 

4. 2  Future  Studies  . 


REFERENCES 


ACKNOWLEDGEMENTS 


I  would  like  to  express  my  appreciation  to  my  committee 
members:  Professors  Julia  N.  Paegle,  Jan  Paegle,  and  Kuo-Nan  Liou 
for  their  support  and  guidance  during  the  preparation  of  this 
thesis.  I  would  especially  like  to  thank  Julia  N.  Paegle  for  the 
many  hours  of  assistance  she  has  given  throughout  my  research. 

I  would  also  like  to  thank  Greg  Dodd,  Jeff  Koenig,  Michael 
Griffen,  Tyree  Wilde,  and  Richard  Petrilla  for  all  the  extra  help 
given  towards  the  writing  of  this  thesis. 

Thanks  should  also  go  out  to  the  rest  of  my  fellow  students 
who  encouraged  me  throughout  my  studies. 

Research  reported  in  this  thesis  has  been  supported  jointly 
by  the  National  Science  Foundation,  the  National  Oceanic  and 
Atmospheric  Administration,  and  the  National  Aeronautics  and  Space 
Administration  under  Grant  Numbers  ATM  8219198  and  NAG  5-127. 


CHAPTER  1 


INTRODUCTION 
1.1  Motivation 

Interest  in  the  global  aspects  of  tropical  meteorology  has 
expanded  markedly  in  the  last  decade.  The  expansion  from  hemi¬ 
spheric  to  global  domains  of  integration  of  numerical  models  in 
operational  meteorological  centers  gave  incentive  to  studies  aimed 
to  improve  our  understanding  of  large  scale  tropical  circulations. 

The  effect  of  tropical  motions  in  extratropical  flows  became  apparent 
as  these  centers  extended  their  forecasts  in  time.  The  impact  of 
tropical  wind  data  in  the  forecasts  of  ultralong  waves  in  the  GLAS/ 
NASA  general  circulation  model  has  been  documented  by  Baker  and 
Paegle  (1983).  They  obtained  large  differences  in  the  tropical  over¬ 
turning  motions  between  forecasts  made  with  and  without  tropical  wind 
data,  suggesting  the  importance  of  heat  driven  circulations  for  global 
scale  predictions.  A  follow  up  study  (Paegle  and  Baker,  1983) 
examined  the  influence  of  tropical  latent  heating  by  comparing  two 
simulations  of  the  GLAS  GCM  started  from  the  same  initial  conditions, 
one  simulation  used  full  physics  and  the  other  contained  no  latent 
heating  in  the  tropical  belt  from  20°  S  to  20°  N.  Their  main 
findings  indicated  the  almost  instantaneous  (within  a  half  day)  impact 
of  tropical  heating  in  the  divergent  tropical  motion  while  the 
rotational  motions  responded  in  about  3  and  5  days  in  tropical  and 


extratropi cal  latitudes. 

These  and  other  studies  have  also  pointed  out  difficulties  in 
properly  simulating  tropical  flows.  This  is  probably  due  to  the 
dominance  of  physical  forcing  in  tropical  latitudes,  in  contrast  with 
mid-latitude  flows  which  are  more  likely  to  be  driven  by  flow  insta¬ 
bilities  and  inertial  effects.  One  primary  source  of  energy  in  the 
tropics  is  the  latent  heat  released  in  precipitation  associated  with 
active  cumulus  clouds  (Holton,  1979). 

Understanding  of  the  physical  processes  which  dominate  in  the 
tropics  requires  quantification  of  different  components  of  the  heat 
budgets  from  complete  data  sets  for  different  seasons  and  time 
periods.  The  current  research  emphasizes  the  latent  heat  component 
in  the  tropical  and  subtropical  region  and  its  effect  in  the  large 
scale  overturning.  The  main  goal  is  to  quantify  the  time  scale  of 
these  motions  and  their  principal  features  as  described  by  the 
heating  rates  obtained  by  GLAS/NASA  and  associated  upper  and  lower 
level  flow  components. 


1.2  Background 

The  Global  Weather  Experiment  (GWE,  December  1,  1978,  November 
30,  1979)  provided  unique  data  sets  to  understand  this  and  other 
characteristics  of  the  tropics.  The  most  comprehensive  program  of 
atmospheric  data  acquisition  ever  attempted  was  successfully 
completed  that  year  with  strong  international  cooperation.  Conven¬ 
tional  observing  systems  were  augmented  utilizing  newly  developed 
techniques.  Research  data  sets  were  compiled  and  archieved  after 
careful  checks  to  remove  observational  or  processing  errors.  These 
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are  denoted  as  level  Il-b  or  1 1 - c  data  sets.  The  assimilation  of  such 
massive  amounts  of  information  presented  a  formidable  challenge. 
Objective  analysis  methods  were  used  to  develop  gridded  fields  using 
the  level  II  data.  These  gridded  fields  are  referred  as  level  1 1 1 - b 
sets  and  are  now  available  from  many  operational  centers  and  research 
institutes.  These  fields  are  somewhat  dependent  on  the  assimilation 
system  used  and  it  is  therefore  important  to  check  the  sensitivity 
of  the  results  to  the  analysis  methods. 

For  this  study,  the  precipitation  data  gathered  during  this 
period  (level  1 1 - c )  are  studied  in  conjunction  with  the  GLAS  heating 
fields  and  circulation  data  obtained  from  the  Ill-b  data  sets. 

Partial  comparisons  of  circulation  statistics  are  presented  using  the 
European  Center  for  Medium  Range  Weather  Forecasting  (ECMWF), 
Geophysical  Fluid  Dynamics  Laboratory  (GFDL) ,  and  Goddard  Laboratory 
for  Atmospheric  Sciences  (GLAS)  analyses. 

Daily  values  of  precipitation  data  were  gathered  globally  during 
1979  directed  to  the  second  objective  of  the  First  GARP  (Global 
Atmospheric  Research  Project)  Global  Experiment  (FGGE):  understanding 
of  basic  physical  forces  affecting  climate.  Data  collected  mainly  as 
a  result  of  this  objective  are  referred  to  as  level  II-c  sets.  The 
global  extent  of  this  enhanced  and  checked  precipitation  data  set  is 
unique  to  the  GWE  year.  It  is  therefore  important  not  only  to 
analyze  this  particular  year,  but  also  to  compare  it  with  more 
extensive  climatologies  whenever  possible  to  isolate  areas  of 
averaged  and  anomalous  atmospheric  behavior.  Such  comparisons  are 
also  included  in  the  present  work. 


CHAPTER  2 


PROCEDURES 

2. 1  Data  Acquisition 

The  following  discussion  describes  analysis  methods  and  data 
acquisition.  The  models  used  to  obtain  the  heat  and  circulation 
fields  are  discussed.  Heating  rates,  horizontal  wind  components, 
velocity  potential,  and  stream  function  fields  are  studied  and  dis¬ 
played  in  this  research.  Information  about  other  fields  is  included 
in  the  discussion  when  necessary  to  complete  the  description  of  the 
physical  phenomena  of  interest. 

2.1.1  Precipitation  Analyses 

The  precipitation  data  is  from  the  FGGE  data  observed  from 
December  1978  -  November  1979.  The  unit  is  tenths  of  millimeters/ 

24  hours  and  there  are  two  categories  of  quality  control: 

1)  accurate,  checked,  and  found  correct;  2)  estimated,  checked,  and 
found  correct.  Category  2)  is  used  to  enhance  the  data  set  which 
otherwise  would  be  very  sparse  in  certain  tropical  areas.  A  trace  of 
precipitation  is  equated  to  zero  amount.  Accumulating  precipitation 
is  averaged  over  the  period  and  used  for  each  day  of  the  accumulating 
period. 

Precipitation  data  are  only  available  over  land.  This  causes 
problems  with  analyses,  contouring,  and  interpretation  of  global 


tropical  data,  due  to  large  oceanic  areas  void  of  data.  Terrain  fea¬ 
tures,  such  as  mountain  versus  valley  observation  points,  island  loca 
tion  and  topographic  features,  observation  density,  and  other  factors 
affect  the  reliability  of  the  data.  Future  research  needs  to  involve 
other  available  information  such  as  satellite  data,  climatology,  etc. 
to  fill  in  data-sparse  regions  and  increase  reliability. 

Recent  developments  in  the  use  and  interpretation  of  polar 
orbiter  satellite  data  makes  possible  obtaining  global  representa¬ 
tion  of  cloud  cover,  cloud  top  temperature,  and  pressure  as  well  as 
surface  temperature  and  vertical  temperature  and  humidity  profiles. 
The  latter  are  routinely  obtained  by  NOAA/NESDIS  for  operational 
purposes  using  the  TIROS  Operational  Vertical  Sounder  (TOVS)  system. 
Novel  methods  of  combining  radiance  data  with  surface  and  atmospheric 
parameters  (Susskind,  et  al.,  1984)  show  promise  for  more  accurate 
global  determination  of  satellite-derived  quantities  and  represent  an 
alternative  to  methods  based  on  statistical  regression  relationships. 
Another  development  measures  water  vapor  over  oceans.  Microwave 
measurements  of  precipitating  versus  nonprecipitating  clouds  were 
derived  by  Hwang,  et  al.  (1984)  using  the  NIMBUS  7  Scanning 
Multichannel  Microwave  Radiometer  (SMMR)  data  from  December  1978 
through  November  1979.  Satellite-derived  precipitation  values 
require  further  study  before  being  used  operationally.  To  be  bene¬ 
ficial,  procedures  that  obtain  surface  rainfall  rates  must  first  be 
developed. 

Precipitation  is  highly  variable  depending  on  local  terrain 
configuration  and  could  attain  very  different  values  at  points 
separated  by  short  distances.  Therefore  only  horizontally  averaged 


values  are  currently  considered. 


2.1.2  Heating  Rate  Fields 

Estimated  of  diabatic  heating  fields  can  be  obtained  in  one  of 
three  ways: 

i)  By  direct  evaluation  of  each  component  of  the  heating 
budget  based  on  parameteri zation  assumptions  (e.g.,  Newell,  et  al., 
1973).  This  method  is  highly  dependent  on  the  parameterization  model 
used. 

ii)  As  residuals  of  the  thermodynamic  energy  equation.  In 
this  approach  individual  terms  are  evaluated  from  observed  data  (e.g. 
Luo  and  Yanai ,  1984)  or  from  gridded  data  (e.g.,  Murakami  and  Ding, 
1982;  Kasahara  and  Mizzi,  1983  using  isobaric  coordinates  and  Wei,  et 
al.,  1983  and  Johnson,  1980  using  isentropic  coordinates). 

iii)  As  residuals  of  the  thermodynamic  energy  equation  where 
individual  terms  are  computed  from  either  the  first  guess  of 
objective  analysis  systems  (e.g.,  Kalnay  and  Baker,  1984)  or  as  fore¬ 
casts  in  numerical  integrations.  In  the  case  of  assimilation  systems 
which  considerably  weaken  the  tropical  overturnings  in  the  course  of 
their  initialization  procedure,  this  method  has  the  clear  advantage 
of  allowing  the  atmosphere  to  spin  up  in  response  to  the  model 
forci ng. 

A  comparison  of  results  obtained  using  method  ii)  by  Kasahara 
and  Mizzi  (1983)  and  the  GLAS  heating  rates  obtained  from  the  first 
guess  analysis  of  their  assimilation  model  (method  iii))  was  done 
Paegle  and  Paegle  (1984). 

Kasahara  and  Mizzi  (1983)  obtained  the  field  as  a  residual  of 
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the  thermodynamic  equations.  The  equation  is: 

?T/3t  +  V-vT  -  wr  =  Q/Cp  (2.1) 

where  r  =  -(T/9)/9g/3p  and  Q/Cp  =  the  residual  heating/cooling  rate. 
The  V  and  T  fields  were  analyzed  every  6  or  12  hours.  The  resolution 
of  this  method  was  1.875  x  1.875  degrees.  The  period  of  study  was 
from  January  27  -  February  11,  1979.  Its  large  scale  features  looked 
realistic,  but  they  isolated  two  possible  shortcomings:  1)  A  semi¬ 
diurnal  wave  #2  oscillation  related  to  the  divergence  field  was  noted 
which  may  be  caused  by  diurnal  variation  of  cumulus  convection,  but 
was  most  probably  due  to  different  data  density  at  06  GMT;  2)  Also, 
an  over  response  in  the  mountainous  areas  was  noted  due  to  model 
divergence  and  vertical  motion  inadequacies  in  the  vicinity  of 
mountains. 

The  GLAS  heating  rates  are  evaluated  in  the  course  of  the  data 
assimilation  using  model  diagnostics  (method  i i i ) ) .  A  damping  time 
step  is  used  for  initialization  to  damp  gt3vity  waves,  instead  of 
the  nonlinear  normal  mode  initialization  used  by  most  operational 
centers.  The  model  resolution  is  4  degrees  x  5  degrees  and  it  has 
full  physics.  The  thermodynamic  equation  used  is: 

t!(Jan  <2'2> 

where  m  =  model  values  and  Q  =  the  difference  between  the  first 

a  n 

guess  and  the  analysis  field.  If  the  first  guess  was  a  perfect  fore¬ 
cast,  fiQan  =  0  and  Qm  =  the  real  atmospheric  heating  rates.  Compari¬ 
son  of  Qm  with  those  obtained  by  Kasahara  and  Mizzi  (1983)  for  the 
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period  January  27  -  February  11,  1979  (Paegle  and  Paegle,  1984) 
showed  that  the  heating  rates  from  the  model  output:  1)  did  not  have 
the  unrealistic  diurnal  oscillations,  and  2)  were  not  as  sensitive 
to  the  divergence  analysis  near  the  mountains  as  method  ii).  The 
GLAS  heating  fields  appeared  excessive  in  certain  locations,  though 
some  computed  peak  heating  rates  of  10°C/day  found  for  this  period 
may  be  realistic.  These  values  are  likely  to  be  observed  in  the 
tropics  where  the  average  value  of  precipitation  is  2  cm  of  precipi¬ 
tation  releasing  approximately  5°C/day  of  heat  over  a  column  of 
atmosphere  (Holton,  1979). 


2.1.3  Circulation  Diagnostics 


The  circulation  variables  have  been  computed  from  three  different 


data  sets,  obtained  by  ECMWF,  GFDL,  and  GLAS/NASA.  Further  comparison 
of  results  obtained  from  different  assimilation  methods  is  desirable 


to  quantify  the  sensitivity  of  conclusions  to  the  analysis  scheme. 

No  heating  rates  were  available  from  the  GFDL  and  ECMWF  assimilation 
cycles,  so  the  intercomparisons  are  limited  to:  U  and  V  wind  compo¬ 
nents,  velocity  potential,  and  stream  function  fields. 


The  GLAS  FGGE  Ill-b  data  set  used  a  successive  correction  method 


which  included  weights  depending  on  the  data  quality  and  density  for 
the  analysis  scheme  (Baker,  1983).  The  analysis  and  forecast  system 
used  the  GLAS  4th  order  GCM  (Kalnay  and  Baker,  1984)  with  a  damping 
time  step  during  initialization  as  previously  indicated  and  included 


all  diabatic  processes.  The  available  analyzed  fields  included 


horizontal  wind  components,  temperature,  humidity,  and  sea  level 


pressure.  Also,  sensible  and  latent  heat,  precipitation,  and  total 


diabatic  heat  fields  were  produced  and  archived. 

The  ECMWF  III -b  data  set  used  intermittent  data  assimilation 
and  a  nonlinear  normal  mode  initialization  (Bengtsson,  1983).  The 
analysis  was  performed  in  pressure  coordinates  while  the  forecast  was 
carried  out  in  a  high  resolution  sigma  coordinate  system.  Its 
analyzed  fields  of  geopotential  heights,  sea  level  pressure,  and 
horizontal  wind  fields  were  not  initialized;  whereas,  the  derived 
fields  of  temperature,  relative  humidities,  and  vertical  velocities 
were  initialized.  The  data  are  available  in  isobaric  surfaces  with 
values  interpolated  under  the  surface  in  mountainous  areas.  The 
vertical  velocities  are  too  weak  as  a  result  of  the  nonlinear  normal 
mode  initialization  (Carr  and  Ramamurthy,  1983;  and  Lambert,  1983). 
These  resulted  in  reduction  of  the  mean  meridional  circulation  and 
wind  divergence  in  the  tropics.  Another  possible  problem  was  the 
pronounced  diurnal  cycle  of  the  tropical  divergence  field.  In  fact, 
the  06  GMT  divergence  was  nearly  one  half  the  values  at  00  GMT  and 
12  GMT  in  certain  regions.  This  may  be  partly  due  to  data  gaps  in 
the  cloud  track  winds  at  06  GMT  which  may  be  obscuring  an  actual 
atmospheric  response  (Paegle,  1983;  and  Julian,  1983). 

All  models  used  the  6  hour  forecast  fields  as  the  first  guess  to 
analyze  new  data.  The  analyses  should  mostly  reflect  observations  in 
regions  of  high  data  density  and  model  forecasts  in  data  sparse 
regions  where  more  of  the  model's  first  guess  was  used.  The 
Geophysical  Fluid  Dynamics  Laboratory  (GFDL)  Level  Ill-b  data  sets 
used  a  continuous  dynamic  assimilation  system  and  a  nonlinear  normal 
mode  initialization.  The  model  was  a  global  spectral  GCM.  The  first 
guess  was  the  most  recent  analysis,  with  weights  assigned  for 


accuracy  and  interpolation  to  fill  in  gaps.  Also,  the  model  results 
versus  analysis  were  weighted  depending  on  reliability  of  the  data. 
As  time  progressed  real  observed  data  were  inserted  into  the  model, 
so  that  the  data  and  the  model's  dynamics  and  physics  were  con¬ 
sistent.  Virtues  of  the  model  were  consistency  among  variables  and 
the  control  of  error  growth.  Problems  were  data  rejection  and  the 
noisy  appearance  of  the  model  (Stern  and  Ploshay,  1983).  Also,  the 
model  underestimated  severe  storm  patterns  at  times  and  its 
vertical  velocities  were  exceedingly  noisy  in  the  tropics  (Miyokada, 
et  al . ,  1983) . 

When  comparing  the  ECMWF  products  to  the  GFDL  model,  the  GFDL 
products  had  shallower  troughs,  stronger  divergent  flows,  and 
stronger  return  flows  in  the  Hadley  Circulation  (S’. urn  and  Ploshay, 
1983). 


2.2  Analysis  Methods 

The  following  discussion  focuses  on  the  format  of  the  precipita 
tion,  heat,  and  circulation  fields,  along  with  some  problems 
involving  these  quantities. 

2.2.1  Precipitation  Analysis 

Precipitation  data  are  only  available  over  land.  Observed 
values  may  vary  greatly  within  short  distances  due  to  regional 
effects.  Therefore  only  horizontally  averaged  data  over  either 
1  7/8 :  x  1  7/8°  grid  size  or  4°  x  5°  grid  are  used  in  this  study. 
These  grid  intervals  are  chosen  to  agree  with  those  of  GFDL,  ECMWF, 


and  GLAS  analyses. 

The  analysis  covers  different  averaging  periods.  Periods 


covered  are  seasonal,  weekly,  and  daily. 

Typical  station  distributions  and  data  coverage  for  the  precipi¬ 
tation  fields  are  shown  in  Figures  2.1  -2.5.  Maximum  number  of 
stations,  using  Figure  2.4,  is  found  in  the  densely  populated  areas  of 
Australia  and  Southeast  Asia.  In  those  regions  the  observation 
counts  are  1792,  420,  and  224  for  grid  areas  centered  at  eastcentral 
Australia,  Malaysia,  and  the  Philippines,  respectively.  The  Hawaiian 
Islands  contribute  about  350  observations  per  grid  point  while  at  the 
maritime  region  of  Indonesia  about  10  observations  per  grid  point 
are  found. 

Figures  2.6  and  2.7  display  averaged  precipitation  amounts. 

Each  grid  point  value  includes  all  reporting  stations  that  have 
acceptable  quality  control  values.  A  look  at  the  seasonal  and  weekly 
graphs  shows  that  the  4°  x  5°  grid  loses  a  lot  of  detail.  It  also 
reduces  the  maximum  amounts  by  averaging  the  precipitation  over  a 
larger  area.  Largest  precipitation  rates  of  3.5  cm/day  are  found 
during  the  summer  season  in  conjunction  with  the  Asian  monsoon  in  the 
higher  resolution  grid.  This  maximum  is  lowered  to  2.2  cm/day  in  the 
4°  x  5°  grid.  The  week  of  January  5-11  displays  a  remarkably  active 
South  Pacific  Convergence  Zone  with  precipitation  amounts  of  about 
6  cm/day  over  the  Pacific  Ocean  and  5.2  cm/day  over  northeastern 
Australia.  These  patterns  are  greatly  distorted  in  the  4°  x  5°  grid 
(Figure  2.5),  with  only  the  Australian  maximum  apparent  in  the 
analysis.  It  is  therefore  preferable  to  use  the  1  7/8°  x  1  7/8°  grid 
to  display  the  precipitation  amounts.  This  is  done  in  the  next 
chapter. 

The  station  density  for  each  grid  varies  as  seen  in  the 
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station  count  figures  mentioned  above.  The  computer  programs  read 
the  data,  perform  any  needed  calculations,  average  values  for  each 
grid  point,  and  arrange  the  individual  grid  values  into  a  world  map 


2.2.2  Heatinq  Rate  Analysis 


Heating  rates  were  derived  by  GLAS/NASA  using  the  methods 
described  in  the  Data  Acquisition  section  presented  above.  Grid 
point  values  are  available  in  the  GLAS  4°  x  5°  analysis  grid  at  9 
sigma  levels.  This  field  includes  latent  heating,  sensible  heating, 
solar  radiation,  and  infrared  cooling.  The  latent  heat  release  is 
of  particular  interest  to  this  study  for  comparison  with  precipita¬ 
tion  data. 

The  heating  fields  obtained  from  GLAS  are  vertically  integrated 
and  averaged  over  different  periods.  Discussion  will  center  on 
seasonal,  weekly,  and  daily  values. 

The  programs  read  the  values  derived  by  the  GLAS  model,  perform 
vertical  averages,  and  arrange  the  individual  grid  values  into  a 
world  map  similar  to  that  obtained  for  the  precipitation  fields. 


2.2.3  Circulation  Analysis 

FGGE  level  Ill-b  data  sets  from  ECMWF  and  GFDL  are  available  in 


1  7/8°  by  1  7/8°  (latitude,  longitude)  grid  as  mentioned  above.  U 
and  V  gridded  fields  are  interpolated  to  a  4°  by  5°  (latitude,  longi¬ 
tude)  grid  for  comparison  with  the  GLAS  data  sets.  A  two  dimensional 
linear  interpolation  scheme  is  used  for  this  purpose. 

Vorticity  (§)  and  divergence  (i)  are  computed  for  each  data  set 


at  several  vertical  levels.  The  Poisson  equations  used  are: 


They  are  numerically  solved  over  the  whole  globe  to  obtain  the 
stream  function  (y)  and  velocity  potential  (v).  The  method  proposed 
by  Paegle  and  Tomlinson  (1975)  is  used  for  this  purpose. 

The  following  chapter  presents  circulation  statistics  for  the 


200  and  850  mb  levels. 


CHAPTER  3 


RESULTS 

3.1  Introduction 

The  main  features  of  the  global  tropical  circulations  may  be 
inferred  from  available  climatologies  (e.g.,  Newell,  et  al.,  1973). 
Tropical  latitudes  provide  a  large  amount  of  the  energy  needed  to 
maintain  atmospheric  circulations  against  frictional  dissipation. 
Sensible  and  latent  heat  are  transported  at  lower  levels  equator- 
ward  and  are  lifted  and  condensed  in  preferred  tropical  regions, 
producing  potential  energy  and  sensible  heat.  Upper  level  circula¬ 
tions  transport  this  energy  poleward.  Condensation  processes  are  of 
great  importance  in  this  process  and  are  the  central  theme  of  the 
current  investigation. 

It  is  not  a  simple  task  to  infer  the  characteristics  of  over¬ 
turnings  from  conventional  data.  Since  precipitation  data  present  a 
number  of  limitations  already  discussed,  it  is  desirable  to  obtain 
heating  rate  estimates  from  other  methods.  The  current  research 
uses  the  heating  rates  derived  at  GLAS/NASA  and  compares  them  with 
precipitation  data  for  cross  validation.  These  comparisons  are  done 
for  a  variety  of  time  and  horizontal  scales:  1)  Seasonal  comparisons 
between  FGGE  data  and  climatological  charts  for  precipitation, 
heating  rates,  and  U  and  V  wind  components  are  presented;  2)  Weekly 
averages  are  contrasted  for  the  two  special  observing  periods  of 


January  5  -  March  1,  1979  and  May  12  -  August  3,  1979;  3)  Heating 

and  velocity  potential  data  are  presented  at  00,  06,  12,  and  18  GMT 
for  the  week  of  January  5  -  11,  1979;  4)  Daily  variations  of  pre¬ 

cipitation  are  compared  with  mid  atmosphere  heating  rate  variations 
for  selected  area  averages  from  January  5,  1979  to  February  28,  1979. 

3.2  Seasonal  Averages 

The  FGGE  year  is  likely  to  remain  the  most  complete  atmospheric 
data  sets  for  the  tropics  and  the  Southern  Hemisphere  for  many  years. 
It  is  important  therefore  to  document  the  differences  observed  with 
respect  to  more  extended  climatologies  to  assess  whether  atypical 
atmospheric  behavior  during  1979  is  due  to  i)  the  enhanced  data  sets, 
or  ii)  time  variations  in  short  climatic  scales  (of  the  order  of  .1 
to  10  years).  Once  anomalous  atmospheric  characteristics  are  isola¬ 
ted  extensive  research  is  needed  to  determine  whether  they  are  due 
to  i)  or  ii)  above.  For  example,  Trenberth  and  van  Loon  (1981) 
and  Trenberth  (1984)  have  discussed  the  rather  anomalous  behavior  of 
the  Southern  Hemisphere  during  FGGE.  In  particular  during  the  months 
of  April  to  July  1979,  they  found  abnormally  low  sea  level  pressures 
over  the  Southern  Hemisphere  (high  pressure  over  the  Northern  Hemi¬ 
sphere).  On  the  other  hand,  Karoly,  et  al.  (1983)  attribute  the  deep 
values  of  the  circumpolar  Southern  Hemisphere  lows  to  the  extended 
network  of  floating  buoys  deployed  during  the  FGGE.  The  discussion 
in  the  present  work  will  concentrate  on  latitudes  from  45°  orth  to 
45°  south. 


The  seasonal  precipitation  averages  display  a  typical  trend  of 
maximum  precipitation  in  the  Southern  Hemisphere  during  the  December 
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through  February  season  (Figure  3.1)  and  in  the  Northern  Hemisphere 
during  the  June  through  August  season  (Figure  3.2).  The  tran¬ 
sition  periods,  March  -  May  and  September  -  November,  are  shown  in 
Figures  3. 1  and  3.2. 

The  FGGE  values  are  compared  with  the  climatology  presented 
in  volume  one  and  two  of  the  book  by  Newell,  et  al.  (1973).  In 
general,  the  dry  and  wet  areas  tend  to  match  up  between  the  FGGE 
data  and  the  climatological  data.  Anomalously  high  precipitation 
values  are  found  over  Morocco,  Cape  York  (Australia),  and 
Madagascar  during  the  winter  season.  During  the  summer,  the  maximum 
precipitation  values  are  almost  twice  as  strong  as  those  found  in 
the  extended  climatology. 

The  winter  season  presents  maximum  values  over  Madagascar, 
Indonesia,  over  the  northeastern  Australian  penninsula  of  Cape  York, 
the  South  Pacific  Convergence  Zone  (SPCZ),  and  the  Amazon  basin.  A 
northward  shift  of  maximum  precipitation  values  is  clearly  apparent 
as  the  year  progresses.  During  the  summer  season  the  strongest 
precipitation  values  are  associated  with  the  Asian  monsoon,  with  the 
maximum  precipitation  rates  located  over  Central  America  in  the 
Western  Hemisphere. 

The  winter  seasonal  precipitation  is  now  compared  to  the  winter 
heating  field  displayed  in  Figure  3.3.  There  is  a  general  corres¬ 
pondence  of  higher  rainfall  with  higher  heating  rates  over  the 
Amazon  basin  and  Indonesia.  Precipitation  maxima  over  Madagascar, 
the  SPCZ  (at  about  15  degrees  south  and  the  dateline),  and  north¬ 
eastern  Australia  have  no  counterpart  in  the  heating  fields.  A 
strong  heating  maxima  of  12.6°C/day  over  the  equatorial  regions  of 
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Figure  3.3.  Precipitation  and  heating  rates. 

(a)  Averaged  precipitation  amounts,  December  1978  -  February 
1979,  1  7/8°  x  1  7/8°  grid,  contour  interval  of  100.0  x 
10-1  mm/day.  Units  are  tenths  of  mm/day. 

(b)  Column  integrated  heating  rates  averaged  for  5  January  - 
1  March,  1979  ( (00Z  +  12Z)/2),  contour  interval  of  50.0 
x  10'l°C/day.  Units  are  tenths  of  degrees  Centigrade/day 


the  Bay  of  Bengal  is  present  in  a  region  with  no  available  precipi¬ 
tation  data.  Heating  fields  include  sensible  heat,  solar  radiation, 
and  infrared  besides  latent  heating  in  the  analyzed  fields.  The 
heating  rates  look  to  be  too  low  over  Australia  and  Northwest  Africa 
and  too  high  over  the  Amazon  basin  of  South  America.  This  possibly 
can  be  attributed  to  the  land  cover  (Trewartha,  1968).  Over 
Australia  and  Northwest  Africa  the  land  surface  has  more  of  a  desert 
type  vegetation  and  a  greater  degree  of  outgoing  radiation,  whereas 
the  Amazon  has  tropical  vegetation  and  less  radiation  escaping  back 
into  the  atmosphere  (Petterssen,  1969).  Therefore,  comparisons  of  the 
vertical  integrated  heat  fields  with  precipitation  data  must  be  done 
with  caution. 

The  winter  season  U  and  V  (Figure  3.4)  components  are  also  com¬ 
pared  to  the  cl imatologi cal  charts  of  both  components  of  the  wind 
velocity.  Main  differences  between  the  FGGE  year  and  the  more 
extended  climatology  are:  1)  The  absence  of  a  separate  westerly 
maxima  over  the  Arabian  penninsula  and  South  America  during  FGGE; 

2)  A  stronger  jet  east  of  Australia  and  a  well-developed  westerly 
maximum  south  of  Africa  during  FGGE.  Both  features  are  found  about 
15  degrees  south  of  their  climatological  location;  3)  The  meridional 
motions  associated  with  the  convective  tropical  regions  of  the 
Southern  Hemisphere  are  in  general  stronger  during  FGGE,  with  north¬ 
ward  flows  about  50«  higher  and  a  southerly  maximum  displayed  over 
Melanesia  of  about  9  m/s  not  present  in  the  more  extended  clima¬ 
tologies. 

Paegle,  et  al.  (1983)  compared  the  global  wind  distributions  of 
the  FGGE  and  Data  System  Test  (DST)  for  the  summer  of  1976.  Their 
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Figure  3.4.  Averaged  components  of  200  mb  wind. 

5  January  -  1  March,  1979,  contour  interval  of  100.0 
x  10*1  m/sec.  Units  are  10'"'  m/sec. 

(a)  U  component 

(b)  V  component 


results  show  stronger  westerly  jets  over  Australia  and  South  Africa 
for  the  FGGE  than  for  the  DST  year.  Similar  conclusions  are 
reached  when  the  FGGE  year  is  compared  to  the  more  extended  clima¬ 
tology  of  Newell,  et  al .  (1973).  These  stronger  circulations  may  be 
related  to  the  Southern  Hemisphere  anomalies  previously  mentioned. 

3.3  Weekly  Averages 

Considerable  interest  has  been  evoked  during  the  last  decade  on 
low  frequency  oscillations  with  30  to  50  day  periods.  These  have 
been  detected  as  eastward  propagating  disturbances  in  the  tropical 
sea  level  presure  field  (Madden  and  Julian,  1971,  1972)  and  eastward 
propagating  planetary  divergence  waves  (Lorenc,  1984).  Meridionally 
propagating  low  frequency  weather  systems  found  by  Krishnamurti  and 
Subrahmanyam  (1982),  especially  over  the  Asian  monsoon  region,  were 
related  to  cloud  cover  by  Yasunari  (1980,  1981)  and  to  active  and 
inactive  spells  of  the  Asian  monsoon  by  Sikka  and  Gadgil  (1980). 
Krishnamurti,  et  al.  (1984)  noted  a  strong  relatinship  between  the 
passage  of  the  divergence  wave  and  wet  and  dry  spells  of  the  monsoon 
region.  The  dependence  of  these  disturbances  on  tropical  convection 
and  divergence  patterns  makes  their  study  directly  relevant  to  our 
current  research. 

These  phenomena  are  best  studied  by  filtering  atmospheric  data 
to  isolate  the  periods  of  interest  as  done,  for  example,  by 
Krishnamurti  and  Gadgil  (1984).  In  this  section  weekly  averages 
are  presented  which  effectively  remove  disturbances  with  short  (less 
than  7  days)  time  scales. 


Results  are  obtained  for  8  weeks  in  the  winter  and  12  weeks  in 
the  summer,  but  only  a  smaller  subset  of  representative  charts  are 


The  precipitation  and  heating  fields  are  shown  in  Figures  3.5 
through  3.11  for  the  first,  fourth,  and  seventh  week  of  the  first 
Special  Observing  Period  (SOP-1)  and  for  the  second,  fifth,  tenth, 
and  twelfth  week  of  the  Second  Observing  Period  (SOP-2).  These 
weeks  are  selected  to  display  the  various  precipitation  regimes  found 
for  different  weeks.  The  precipitation  and  heating  fields  show  an 
overall  trend  of  increased  precipitation  corresponding  to  areas  of 
higher  heating  rates.  This  is  especially  noticeable  when  winter  and 
summer  weeks  are  contrasted.  As  the  highest  areas  of  precipitation 
move  from  the  Southern  Hemisphere  in  the  winter  to  the  Northern 
Hemisphere  in  the  summer  months,  so  do  the  heating  fields.  This  is 
also  a  response  to  the  incoming  radiation  flux  moving  to  the 
Northern  Hemisphere  as  the  seasons  progress  from  winter  to  summer. 

The  first  week  in  SOP-1  is  characterized  by  a  very  active  South 
Pacific  Convergence  Zone  with  values  up  to  6  cm/day  observed  in 
equatorial  latitudes  close  to  the  dateline.  Three  weeks  later  the 
maximum  heating  rates  are  found  over  the  Indian  Ocean,  rainfall  of 
8  cm/day  are  reported  over  the  northeastern  coast  of  Australia,  while 
precipitation  has  markedly  decreased  over  the  dateline.  Maximum 
values  are  found  again  close  to  the  dateline  3  weeks  later. 

These  features  will  be  related  to  circulation  characteristics  below. 

The  strongest  summer  precipitation  values  are  found  for  the 
twelfth  week  of  SOP-2  with  values  of  8  cm/day  over  India.  The  one 
week  averages  for  May  and  June  (Figures  3.8  and  3.9)  on  the  pre¬ 
onset  phase  of  the  Asian  monsoon,  display  typical  precipitation 
maxima  over  southeast  China,  Burma,  and  the  Bay  of  Bengal.  After 
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Figure  3.5.  Precipitation  and  heating  rates  for  5  January 
11  January,  1979. 

(a)  Averaged  precipitation  amounts,  1  7/8°  x  1  7/8° 
grid,  contour  interval  of  100.0  x  10_1  mm/day. 
Units  are  tenths  of  mm/day. 

(b)  Averaged  column  integrated  heating  rates  ( ( 1 00Z 
+  1 2Z ) / 2 ) ,  contour  interval  of  50.0  x  10_l°C/day 
Units  are  tenths  of  degrees  Centigrade/day. 
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Figure  3.7.  Precipitation  and  heating  rates  for  16  February 
22  February,  1979. 

(a)  Averaged  precipitation  amounts,  1  7/8  x  1  7/8 
grid,  contour  interval  of  100.0  x  10-'  mm/day. 
Units  are  tenths  of  mm/day. 

(b)  Averaged  column  integrated  heating  rates  ( ( 00Z  + 
12Z)/d),  contour  interval  of  50.0  x  10_'°C/day. 
Units  are  tenths  of  degrees  Centigrade/day. 


UtUud*  Latitude 


f 

r. 

f 

f 


If 

>L. 


Figure  3 
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.8.  Precipitation  and  heating  rates  for  12  May  - 
1979. 

Averaged  precipitation  amounts,  1  7/8°  x  1  7/8° 
contour  interval  of  100.0  x  10"^  mm/day.  Units 
tenths  of  mm/day. 

Averaged  column  integrated  heating  rates  ( ( 1 00Z 
12Z/2),  contour  interval  of  50.0  x  l0"1°C/day. 
Units  are  tenths  of  degrees  Centigrade/day. 
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Figure  3.9.  Precipitation  and  heating  rates  for  2  June  - 
8  June,  1979. 

(a)  Averaged  precipitation  amounts,  1  7/8°  x  1  7/8'  grid, 
contour  interval  of  100.0  x  10-1  mm/day.  Units  are 
tenths  of  mm/day. 

(b)  Averaged  column  integrated  heating  rates  ( ( 00Z  + 
12Z)/2),  contour  interval  of  50.0  x  10'1'C/day. 

Units  are  tenths  of  degrees  Centigrade/day. 
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Figure  3.10.  Precipitation  and  heating  rates  for  7  July 
13  July,  1979. 

(a)  Averaged  precipitation  amounts,  1  7/8°  x  1  7/8° 
grid,  contour  interval  of  100.0  x  10"!  mm/day. 
Units  are  tenths  of  mm/day. 

(b)  Averaged  column  integrated  heating  rates  ( ( 00Z 

1 27 ) / 2 ) ,  contour  interval  of  50.0  x  10_'°C/day. 
Units  are  tenths  of  degrees  Centigrade/day. 
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Figure  3.11.  Precipitation  and  heating  rates  for  28  July  - 
3  August,  1979. 

(a)  Averaged  precipitation  amounts  1  7/8°  x  1  7/8°  grid, 
contour  interval  of  100.0  x  10""'  mm/day.  Units  are 
tenths  of  mm/day. 

(b)  Averaged  column  integrated  heating  rates  ((00Z  + 
12Z/2),  contour  interval  of  50.0  x  10_l°C/day. 

Units  are  tenths  of  degrees  Centigrade/day . 


June  18  the  Asian  summer  monsoon  of  1979  was  fully  established. 

Heavy  rainfalls  are  recorded  over  India  from  this  date  on,  with 
intermittent  wet  and  dry  spells.  These  features  are  more  closely 
seen  in  Figures  3.12  through  3.14  which  show  the  week  averages  sub¬ 
tracted  from  the  SOP-1  and  SOP-2  averages. 

High  precipitation  rates  over  Australia  in  the  winter  time  have 
no  counterparts  in  the  heating  rates.  Also,  over  South  America  the 
high  heating  rates  do  not  have  as  high  precipitation  values  as  other 
areas  of  comparable  high  heating  rates.  The  Australian  differences 
could  be  due  to  the  lack  of  thick  vegetation  and  corresponding 
greater  radiational  cooling.  Over  South  America,  high  precipitation 
values  may  be  missed  by  the  sparsity  of  the  reporting  precipitation 
network.  This  is  especially  true  over  vast  ocean  areas  where  at 
times  there  are  not  data.  The  Indian  Ocean  maxima  could  not  be 
validated  with  precipitation  data  due  to  the  lack  of  stations  in  this 
area. 

Velocity  potential  and  stream  function  weekly  averages  as  well  as 
their  differences  from  seasonal  averages  are  shown  in  Figures  3.15  - 
3.20.  These  fields  determine  the  rotational  and  divergent  wind 
components  through  the  relationships : 

V  =  k  X  Vv  (3.1; 

vd1,  ■  •*  •  <3-2> 

The  strength  of  the  wind  field  is  determined  by  the  spacing  of 
isolines  in  both  cases.  Though  global  fields  are  presented  the 
discussion  will  center  in  the  tropical  latitudes  of  interest  to  this 
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Figure  3.12.  Averaged  precipitation  amounts  differenced  from 
period  average  (5  January  -  1  March,  1979),  1  7/8°  x  1  7/8°  grid 
contour  interval  of  100.0  x  10"1  mm/day.  Units  are  tenths  of 
mm/day,  January  to  February,  1979. 

(a)  5  January  -  11  January,  1979 

(b)  26  January  -  1  February,  1979 

(c)  16  February  -  22  February,  1979 
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Figure  3.15.  Averaged  fields,  5  January  -  11  January,  1979 
Parallel  hatched  areas  represent  maximum  nondivergent  flow 
strength.  Arrows  represent  direction  of  divergent  or  non¬ 
divergent  flow. 

(a)  Stream  function,  200  mb,  contour  interval  of  25.0 
x  106  m2/sec.  Units  are  106  m2/sec. 

(b)  Velocity  potential,  200  mb,  contour  interval  of 
25.0  x  105  m2/sec.  Units  are  10^  m2/sec. 

(c)  Velocity  potential,  850  mb,  contour  interval  of 
25.0  x  lO^  m2/sec.  Units  are  10$  m2/sec. 
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Figure  3.16.  Averaged  fields  differenced  from  period 
average  (5  January  -  1  March,  1979),  5  January  -  11 
January,  1979.  Parallel  hatched  areas  represent  maximum 
nondivergent  flow  strength.  Arrows  represent  direction 
of  divergent  and  nondivergent  flow. 

(a)  Stream  function,  200  mb,  contour  interval  of 
50.0  x  10®  m^/sec.  Units  are  10®  m2/sec. 

(b)  Velocity  potential,  200  mb,  contour  interval 
of  25.0  x  10$  m2/sec.  Units  are  10®  m2/sec. 

(c)  Stream  function,  850  mb,  contour  interval  of 
50.0  x  10®  m^/sec.  Units  are  10®  m^/sec. 
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Figure  3.17.  Averaged  fields,  26  January  -  1  February, 
1979.  Parallel  hatched  areas  represent  maximum  non- 
divergent  flow  strength.  Arrows  represent  direction  of 
divergent  or  nondivergent  flow. 

(a)  Stream  function,  200  mb,  contour  interval  of 
25.0  x  106  m^/sec.  Units  are  106  m2/sec. 

(b)  Velocity  potential,  200  mb,  contour  interval 
of  25.0  x  105  m2/sec.  Units  are  10$  m2/sec. 

(c)  Velocity  potential,  850  mb,  contour  interval 
of  25.0  x  10$  m^/sec.  Units  are  105  m^/sec. 
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Figure  3.18.  Averaged  fields  differences  from  period 
average  (5  January  -  1  March,  1979),  26  January  -  1  February, 
1979.  Parallel  hatched  areas  represent  maximum  nondivergent 
flow  strength.  Arrows  represent  direction  of  divergent  and 
nondivergent  flow. 

(a)  Stream  function,  200  mb,  contour  interval  of  50.0 
x  105  m^/sec.  Units  are  10^  m2/sec. 

(b)  Velocity  potential,  200  mb,  contour  interval  of 
25.0  x  105  m2/sec.  Units  are  105  m2/sec. 

(c)  Stream  function,  850  mb,  contour  interval  of 
50.0  x  105  m2/sec.  Units  are  106  m^/sec. 


Figure  3.19.  Averaged  fields,  16  February  -  22  February,  1979 
Parallel  hatched  areas  represent  maximum  nondivergent  flow 
strength.  Arrows  represent  direction  of  divergent  or  non¬ 
divergent  flow. 

(a)  Stream  function,  200  mb,  contour  interval  of 
25.0  x  106  rn^/sec.  Units  are  106  m2/sec. 

(b)  Velocity  potential,  200  mb,  contour  interval  of 
25.0  x  10$  m2/sec.  Units  are  105  m2/sec. 

(c)  Velocity  potential,  850  mb,  contour  interval  of 
25.0  x  1 05  m2/sec.  Units  are  10$  m^/sec. 
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Figure  3.20.  Averaged  fields  differenced  from  period 
average  (5  January  -  1  March,  1979),  16  February  -  22 
February,  1979.  Parallel  hatched  areas  represent  maximum 
nondivergent  flow  strength.  Arrows  represent  direction  of 
divergent  and  nondivergent  flow. 

(a)  Stream  function,  200  mb,  contour  interval  of 
50.0  x  10$  m^/sec.  Units  are  10^  m^/sec. 

(b)  Velocity  potential,  200  mb,  contour  interval  of 
25.0  x  10^  m2/sec.  Units  are  1 05  m2/sec. 

(c)  Stream  function,  850  mb,  contour  interval  of 
50.0  x  105  m2/sec.  Units  are  106  m2/sec. 


study.  Divergence  aloft  and  convergence  in  lower  levels  correspond 
to  higher  tropical  precipitation  and  heating  values.  At  extra- 
tropical  latitudes  the  contribution  from  latent  heat  release  to  the 
total  heating  field  is  small. 

The  week  of  January  5-11  is  characterized  by  strong  divergent 
motions  emanating  from  the  central  Pacific  and  South  Pacific 
Convergence  Zone  where  precipitation  and  heating  maxima  are 
located.  Another  divergence  center  is  located  over  South  America. 
Strong  subtropical  westerly  jets  in  the  Northern  Hemisophere  are 
found  at  similar  longitudes.  Sampson  (1982)  documented  similar 
relationships  in  his  weekly  analysis  of  F6GE  data.  The  200  mb 
anomalies  (defined  as  differences  from  seasonal  averages)  indicate 
that  these  two  divergence  centers  are  stronger  than  the  seasonal 
averages  and  are  found  in  general  in  conjunction  with  equatorial 
westerlies.  The  week  of  January  26  -  February  1  displays  divergent 
circulation  over  South  America  and  the  dateline  close  to  seasonal 
values.  The  main  divergent  flow  anomaly  is  centered  over  the  Indian 
Ocean  in  apparent  response  to  the  heating  anomaly  previously 
discussed.  Equatorial  easterlies  are  found  in  the  200  mb  anomaly 
chart  over  and  to  the  west  of  the  velocity  potential  minimum  with 
westerlies  to  the  east.  The  week  of  February  16-22  is  similar  to 
the  seasonal  averages  over  the  Southern  Hemisphere  tropics,  while  a 
divergence  anomaly  is  apparent  over  Nigeria. 

These  pulsations  of  tropical  convection  and  related  overturnings 
may  be  partly  due  to  a  long  periodicity  perturbation  (such  as  the 
30  -  50  day  wave)  propagating  through  the  global  tropics. 


3.4  Daily  Variations 

Rainfall  rates  are  compared  with  heating  rates,  integrated  over 
a  column  from  about  250  to  600  mb  every  6  hours.  Heat  changes  are 
plotted  for  every  tenth  degree  centigrade  and  rainfall  rates  are 
plotted  for  every  mi llimeter/day.  The  selected  mid  atmospheric 
heating  rates  are  chosen  to  better  represent  latent  heat  release. 

Five  different  areas  are  plotted.  The  five  areas  investigate 
portions  of  South  America  (Figure  3.21),  the  Indian  Ocean  (Figure 
3.21),  the  South  Pacific  (Figure  3.21),  Australia  (Figure  3.22), 
and  Africa  (Figure  3.22).  Area  boundaries  are:  South  America  - 
4  degrees  south,  21  degrees  south.  43  degrees  west,  58  degrees  west; 
the  Indian  Ocean  -  4  degrees  north,  11  degrees  south,  49  degrees 
east,  79  degrees  east;  the  South  Pacific  -  11  degrees  south,  24 
degrees  south,  142  degrees  west,  167  degrees  west;  Australia  -  4 
degrees  south,  17  degrees  south,  131  degrees  east,  152  degrees  east; 
and  Africa  -  4  degrees  south,  21  degrees  south,  23  degrees  east, 

43  degrees  east.  They  are  plotted  over  a  55  day  period. 

Variability  increases  in  one  field  are  related  to  greater  vari¬ 
ability  in  the  other  field.  Both  variables  appear  to  increase  and 
decrease  in  a  similar  fashion.  The  largest  precipitation  rates  are 
found  over  the  SPCZ  during  the  first  2  weeks  of  SOP-1.  This  area 
shows  higher  heating  rates  when  there  are  higher  precipitation 
rates.  Other  areas  do  not  exhibit  the  dramatic  long  term  change  from 
high  to  low  rates  within  the  55  day  period,  as  does  the  South  Pacific 
area,  though  other  oscillations  with  time  scales  of  1  to  2  weeks  can 
be  detected  in  these  graphs.  The  heating  rates  have  a  strong  semi¬ 


diurnal  component. 
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Figure  3.22.  Heating  and  rainfall  rates  versus  time, 

5  January  -  28  February,  1979,  Australia  and  Africa. 

(a)  Australia,  257  to  300  reporting  stations  per  day 

(b)  Africa,  40  to  67  reporting  stations  per  day 


The  differences  between  the  circulation  models  and  the  differ¬ 
ences  from  00  GMT  are  displayed.  The  GLAS  heating  rates  at  00  GMT 
and  the  differences  from  00  GMT  are  shown  on  Figure  3.23  and 
velocity  potentials  from  ECMWF,  GFDL,  and  GLAS/NASA  analysis  models 
are  displayed  in  Figures  3.24  -  3.26. 

The  high  heating  rates  over  South  America  do  not  have  the 
intense  divergence  at  200  mb  as  other  high  heating  rate  areas  pro¬ 
duced  by  these  models.  All  the  models  match  areas  of  upper  level 
divergence  with  higher  heating  rates,  and  areas  of  upper  level  con¬ 
vergence  with  lower  heating  rates. 

The  heat  field  differences  from  00  GMT  show  an  apparent  move¬ 
ment  with  time.  The  positive  heating  rates  follow  the  daylight 
hours.  Land  masses  probably  influence  this  diurnal  trend  by  dis¬ 
torting  heating  rates  toward  positive  values  over  longitudes 
corresponding  to  these  areas.  The  flux  of  incoming  and  outgoing 
radiation  seems  to  assert  a  greater  effect  on  the  heating  rates. 

This  effect  could  also  be  due  to  extended  areas  of  increased  daytime 
cumulus  convective  activity,  especially  over  land  masses,  but  should 
be  more  localized  than  seen  on  the  charts.  Since  the  GLAS  heating 
rate  differences  and  GLAS  velocity  potential  differences  are  from 
the  same  model  they  show  a  relationship;  whereas  the  ~ther  two  models 
do  not.  This  relationship  indicates  the  consistency  between  the  first 
guess  of  the  model  and  the  actual  data.  The  flow  intensity  changes 
of  the  divergence  field  aligns  mostly  along  and  inside  the  border  of 
the  positive  heating  area  (Figures  3.23,  3.24,  and  3.26).  Most  areas 
decrease  the  original  divergent  flow  in  response  to  the  boundaries  of 
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cooling  and  heating.  It  must  be  remembered  there  is  a  correspon¬ 
dence  between  the  precipitation  fields,  the  velocity  potentials 
(divergence),  and  heating  rates,  since  latent  heat  release  is  a 
component  of  the  heating  rates. 

The  individual  velocity  potential  patterns  at  00  GMT  (Figures 
3.24  -  3.26)  of  the  three  models  are  very  similar,  but  when  intensi¬ 
ties  are  viewed,  one  model  is  much  different  than  the  others.  The 
GLAS  and  GFDL  data  look  quite  similar,  but  the  ECMWF  values  are  much 
smaller.  The  center  values  in  the  GLAS  and  GFDL  models  are  quite 
close  to  each  other.  The  GLAS  and  the  ECMWF  fields  agree  where  the 
two  lobes  of  the  velocity  potential  come  off  the  deepest  low  center; 
the  GFDL  field  does  not.  The  GLAS  and  ECMWF  fields  emphasize  the 
northern  lobe  and  GFDL  fields  emphasize  the  western  lobe.  These 
results  agree  with  results  mentioned  in  the  circulation  acquisition 
section.  That  is,  the  GFDL  model  has  stronger  divergent  flow  and 
shallower  troughs  compared  to  the  ECMWF  model. 

The  differences  from  00  GMT  (Figures  3.24  -  3.26)  show  GLAS 
values  having  a  much  greater  change  form  the  00  GMT  data  than  the 
other  two  fields.  The  ECMWF  data  have  the  smoothest  appearance.  The 
change  of  the  divergent  flow  with  time  is  another  apsect  of  these 
differences.  The  GLAS  model  has  much  greater  flow  intensity  changes 
of  velocity  potential  than  the  other  two  models  (Figures  3.24  -  3.26), 
with  intensity  changes  occurring  in  different  areas.  Also,  the 
continents  can  influence  the  difference  fields.  Good  examples  are 
the  GLAS  differenced  fields  (Figure  3.26)  over  South  America  at 
06  GMT,  Africa  at  12  GMT,  and  the  east  coast  of  North  America  at 
18  GMT.  The  other  models  show  the  same  relationship  having  most 


areas  either  over  continents  or  divided  by  their  coastlines.  Also, 
there  is  a  diurnal  trend  towards  a  decrease  or  increase  of  divergent 
flow  favoring  oncontinental  200  mb  flow  late  evening  to  early 
morning  and  offcontinental  200  mb  flow  during  the  afternoon  and  early 
evening,  for  example,  over  Asia  and  South  America.  Between  those 
times  the  changes  in  divergent  flow  and  area  of  change  decrease. 

Here,  again,  all  three  fields  have  similar  appearances.  Overall, 
the  GLAS  field  has  greater  change  with  time,  the  ECMWF  data  have  much 
lower  intensities,  and  the  GFDL  field  has  shallower  troughs. 


CHAPTER  4 


CONCLUSIONS 

4.1  Data 

The  purpose  of  this  study  is  to  take  the  FGGE  precipitation 
data  and  to  put  it  into  a  usable  form.  This  is  accomplished  by 
using  data  that  are  without  error,  averaging  that  data  over  time  and 
space,  and  arranging  the  grids  of  data  on  a  world  format.  This  is 
done  from  45°  north  to  45°  south  for  daily  averages  over  various 
periods  (seasonal,  weekly,  and  daily). 

The  data  are  compared  to  climatology  and  other  fields 
created  from  the  FGGE  data.  Although  it  has  been  noted  that  the  GWE 
year  was  atypical  (Paegle,  1984),  the  results  here  show  a  general 
trend  of  correspondence  between  the  data.  Of  course,  there  are 
regions  of  disagreement  and  uncertainty.  Some  of  the  uncertainty  can 
come  from  the  precipitation  data  being  quite  different  from  precipi¬ 
tation  values  a  short  distance  away  but,  overall,  the  values  are 
fairly  representati ve,  especially  over  areas  of  dense  observations. 
Another  direct  problem  is  the  scarcity  of,  or  nonexistence  of,  data 
in  certain  locations,  especially  over  the  oceans. 

The  seasonal  charts  show  northward  movement  of  rainfall  maxima 
from  winter  to  summer,  with  higher  than  normal  values  located  over 
Morocco,  Cape  York  (Australia),  and  Madagascar.  Higher  heating  rates 
correspond  to  higher  rainfall  rates,  although  heating  rates  are  low 


over  Madagascar,  SPCZ,  northwest  Africa,  and  northeast  Australia,  and 
high  over  the  Amazon  basin.  The  U  and  V  wind  components  have  no 
western  maxima  over  the  Arabian  penninsula  and  South  America,  have 
a  stronger  jet  east  of  Australia  and  south  of  Africa,  both  15  degrees 
further  south,  and  have  stronger  meridional  motions  associated  with 
convective  tropical  regions,  when  compared  to  extended  climatology. 

The  weekly  charts  have  an  overall  trend  of  high  rainfall  rates 
over  high  heating  rates,  with  the  same  correspondence  as  the 
seasons.  Values  of  up  to  8  cm/day  are  recorded  over  India  and  the 
northeast  coast  of  Australia.  Also  identified  is  a  3  week 
oscillation  between  high  rainfall  rates  over  the  Indian  Ocean  and 
dateline.  The  divergence  aloft  corresponds  with  higher  heating  rates 
and  rainfall  rates,  and  oscillates,  as  above. 

The  daily  variations  for  selected  areas  show  heating  rates 
having  semidiurnal  fluctuations.  As  the  variability  of  the  heating 
rates  increase,  the  rainfall  rates  increase.  These  variables  appear 
to  increase  and  decrease  together. 

The  hourly  heating  rates,  rainfall  rates,  and  convergence  again 
correspond  as  above  relations.  There  is  an  apparent  movement  of 
positive  heating  values  with  daylight  hours,  although  distorted  over 
land  masses.  Divergent  flow  changes  occur  mainly  along  and  inside 
the  positive  value  area  in  the  GLAS  model.  The  increase  or  decrease 
of  the  velocity  potential  flow  also  changes  diurnally.  The  daylight 
hours  favor  offcontinental  200  mb  flow,  whereas  the  late  evening  to 
early  morning  hours  favor  oncontinental  200  mb  flow.  This  also  is  an 
apparent  influence  of  the  velocity  potential  differences  by  the  con¬ 
tinents.  The  ECMWF  model  has  lower  divergence  intensities,  the  GLAS 


model  has  greater  change  with  time,  and  the  GFDL  model  has  shallower 
troughs. 

Some  of  the  regional  differences  are  also  probably  attributable 
to  the  following  factors:  land  cover  (which  affects  the  amount  of 
radiation  retained),  latitude  (which  determines  which  heat  component 
is  dominant),  amount  of  nonprecipitating  clouds,  observational 
density,  sea  surface  temperatures ,  and  many  others.  These  regional 
differences  need  to  be  studied  in  more  depth  in  the  future. 

4.2  Future  Studies 

More  can  be  done  to  this  material  to  make  it  more  complete.  Th 
use  of  satellite  data  is  easily  accessible.  The  TIROS  Operational 
Vertical  Sounder  (TOVS)  system  can  be  used  to  get  global  coverage 
(e.g.,  Susskind,  et  al.,  1984).  A  recent  development  in  using  the 
microwave  measurement  of  water  vapor  is  yielding  very  good  results 
over  the  oceans.  It  can  measure  precipitable  water  versus  nonpre¬ 
cipitating  clouds.  For  these  results  the  NIMBUS  7  Scanning  Multi¬ 
channel  Microwave  Radiometer  (SMMR)  data  derived  by  Hwang, 

et  al.  (1984)  could  be  used  for  precipitation  data  over  the 
oceans  if  future  research  developments  provide  surface  rainfall 
rates.  Due  to  the  sparseness  of  the  data  over  the  oceans,  the 
increase  of  this  information  would  greatly  add  to  the  data  input 
from  these  areas  on  future  studies.  The  data  were  also  gathered  for 
the  period  from  December  1978  through  November  1979,  the  same  time 
frame  as  the  FGGE  period. 
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